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^ ! Abstract 

g I We perform the Monte Carlo simulations of the hard-sphere lattice gas on the 

I ' simple cubic lattice with nearest neighbour exclusion. The critical activity 

pH ■ is estimated, Zc = 1.0588 ± 0.0003. Using a relation between the hard-sphere 

O ■ lattice gas and the antiferromagnetic Ising model in an external magnetic 

field, we conclude that there is no re-entrant phase transition of the latter 
on the simple cubic lattice. 



1 Introduction 



The antiferromagnetic Ising model shows a phase transition in an external 
magnetic field but the ferromagnetic one has a critical point only in zero 
field. The critical line surrounds the antiferromagnetic ordered phase ||l|, §]. 
A Hamiltonian is 

{ij) « 

where Sj is an Ising spin variable located zth lattice site and which takes 
on the value +1 and —1. The first summation is over all nearest neighbour 
pairs on a lattice and the second over all lattice sites. J(< 0) is the exchange 
interaction. H is an external magnetic field. Many authors have studied the 
system by various methods: Bethe approximation 0], mean field approxi- 
mation 0], constant coupling approximation Kikuchi approximation [^], 
series expansions |[ , Monte Carlo simulations [1^, 0, |12], |T3|, [l^ , trans- 
fer matrix [Q, renormalization group with periodic cell clusters |ll5| , |16 
phenomenological renormalization group and transfer matrix W7\, [T3], exact 
calculations of an interfacial model [jl9|, finite-size scaling and transfer matrix 



[ Pm , analyses of a 16- vertex model PT |. 

On a bipertite lattice there is a critical field, He = q\J\ where g is a coor- 
dination number. At T = all of the nearest neighbour pairs is antiparallel 
for \H\ < He and is parallel for \H\ > He- At if = the antiferromagnetic 
Ising model is equivalent to the ferromagnetic one. Then in the former a 
phase transition occurs at the Neel temperature, Tn, which is equal to the 
critical temperature of the latter. The three points, (T, H) = (0, He), (Tn, 0), 
and (0, —He), are on the critical line. Since the phase diagram is symmetrical 
for H = 0, we consider non-negative H case. 

When the slope of the critical line at (T, H) = (0, He) is positive, a re- 
entrant phase transition occurs. The system is in the paramagnetic, the 
antiferromagnetic ordered, and the paramagnetic phase as T is decreased 
when H is slightly above He- It is negative on the square lattice 
and is positive on the body-centred cubic lattice 0, |12|, |18 
simple cubic lattice there is no conclusive result |I8|. The purpose of this 
paper is to identify whether there is the re-entrant phase transition of the 
antiferromagnetic Ising model on the simple cubic lattice or not. 

We carry out Monte Carlo simulations of a hard-sphere lattice gas rather 
than the antiferromagnetic Ising model. In the next section we review a 
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relation between them. We describe a Monte Carlo algorithm for the hard- 
sphere lattice gas and define physical quantities measured in section |^. In 
section |^ we present Monte Carlo results. A summary is given in section ^. 



2 Hard-sphere lattice gas 



The slope of the critical line relate to the critical activity of a hard-sphere 
lattice gas with nearest neighbour exclusion. The atoms can occupy lattice 
sites under two conditions: only one atom may occupy a lattice site and 
they interact with infinite repulsion of nearest neighbour pairs. There are 
many studies of the system: series expansions [E^, El, El, E3, 12H], finite-size 



scaling and transfer matrix ||2^, Bethe and ring approximations [|28[, transfer 
matrix [H, H, H, corner transfer matrix and series expansions ||32|, exact 
calculations 



, Monte Carlo simulations 



The grand partition function is 



N 



N 



ZviN), 



(1) 



where z is an activity and V is the number of lattice sites. Zv{N) is the 
number of configurations in which there are N atoms in a lattice of V sites. At 
z = +00 a ground state configuration is that the atoms occupy all the sites of 
one sublattice and the other is vacant. There is no atom at z = 0. The critical 
activity, Zc, relates to the slope, a*, of the critical line at {T,H) = {0,Hc) 
defined hj H = + a^fegT [0, 



2|, H 



Inzc 



(2) 



A result by series expansions is that Zc = 1.09(7) for the simple cubic lattice 
2^. It shows that the value of a* is negative with a large error; a* = —0.04(3) 



IS 



As far as we know, a Monte Carlo study of a hard-sphere lattice gas has 
been published only on the square lattice |Q. For the first time we perform 
Monte Carlo simulations on the simple cubic lattice. If we use a Monte 
Carlo method for the antiferromagnetic Ising model, we must simulate the 
system at very low temperature to estimate the slope of the critical line at 
(T, H) = (0, He). It is difficult to get rehable data. 
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3 Monte Carlo simulations 



We use the Metropolis Monte Carlo technique ^ to simulate the hard- 
sphere lattice gas (|l]) on the simple cubic lattice of V sites, where V = 
L X L X L (L = 10, 12, 14, 16, 18, 20, 22, and 24), under fully periodic 
boundary conditions. 

According to Meirovitch [^], we adopt the grand canonical ensemble. 
The algorithm is as follows. 

1. Choose a lattice site. 

2. If its nearest neighbour sites are occupied, go to |l|. If they are vacant, 
go to the next step. 

3. Generate a random number, r G [0, 1]. 

4. When an atom occupy the lattice site, we remove it if r < min(z~^, 1). 
When it is vacant, we add an atom if r < mm{z, 1). 

5. Goto|. 

We start each simulation from a large activity with a ground state con- 
figuration and then gradually decrease an activity. The pseudorandom num- 



bers are generated by the Tausworthe method We measure physical 

quantities at an activity over 10^ Monte Carlo steps per site after discarding 
5 X 10"^ Monte Carlo steps per site to attain equilibrium. We have checked 
that simulations from the ground state configuration and the configuration 
without atoms gave consistent results and there was no hysteresis. Each run 
is divided into ten blocks. Let us the average of a physical quantity, O, in 
each block (O)^; i = 1, 2, . . . , 10. The expectation value is 

1 1° 



i=l 

The standard deviation is 
Let us define a density by 

p = N/V 
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where N is the number of the atoms in the lattice of V sites and an order 
parameter by 

R = 2{Na-Nb)/V 

where A'a (^b) is the number of the atoms in the A (B)-sublattice and 
= A^A + Ab- We measure the isothermal compressibility: 



K = v{{p^)-{py)/{p), (3) 

the staggered compressibility: 

= V{{R^) - {\R\y)/4, (4) 
and the fourth-order cumulant of R p9| : 

U = 1-^W)/{W- (5) 

4 Monte Carlo results 

Figure 1 shows the activity dependence of the isothermal compressibility, 
defined by (H) for various lattice sizes. The solid curves are obtained 



by the smoothing procedure of the i?-spline \W]. As L increases, the shape 
of the curve becomes sharper. In figure 2 we show the activitiy dependence 
of the staggered compressibility, x^Li^)y defined by (^) for various lattice 
sizes. It does not seem that the position of the peak shifts to contrast those 
of kl{z). We show the activity dependence of the fourth-order cumulant, 
Ul{z), of R defined by (|]) for various lattice sizes in figure 3. There is an 
intersection between the curves with the size L and L + 2. The positions of 
these intersections are within a narrow region. 

We define effective critical activities, z^^^{L) and z^^^{L), as the peak 
position of kl{z) and xll^); respectively, and ^^^^^(L) by 

They will converge to the critical activity, ^c, as L — > +oo. We plot them 
against 1/L in figure 4. 

We decide to estimate Zc from z'^^J^L) by the following reasons. Al- 
though z^g^^(L) seems to behave systematically for L, it is difficult to ex- 
trapolate Zc from it since we need a precise value of a critical exponent v: 
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^maxW ~ -2^c ~ L^^^" . We cannot see systematic behaviour for L in z'^^-^{L) 
and 2;^ogg(L). The latter is more scattering than the former. We get the 
result, Zc = 1.0588(3), by the arithmetic mean from the data except those of 
L = 10. Using (0), we have that a* = —0.02857(14). These are consistent 
with previous results, = 1.09(7) p6[| and a* = —0.04(3) |]l8l, within errors. 
Our results are more accurate than theirs. 



5 Summary 

We performed the Monte Carlo simulations of the hard-sphere lattice gas 
on the simple cubic lattice. We estimated the critical activity, Zc, to be 
1.0588 ± 0.0003. It means that the slope of the critical line at {T,H) = 
(0, He) is negative in the antiferromagnetic Ising model on the simple cubic 
lattice; a* = —0.02857 ± 0.00014. We conclude that there is no re-entrant 
phase transition in this system. 

In closing this paper we want to mention that the antiferromagnetic Ising 



model is equivalent to a spin-one Ising model Our result confirms that 
a re-entrant phase transition occurs in the latter on the simple cubic lattice. 

We have carried out the simulations on two personal computers with the 
486DX2/66MHZ CPU and the Linux operating system (SLS 1.0.3 + JE 0.9.3, 
Slackware 2.0.0 + JE 0.9.5). 
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Figure captions 



Figure 1 Activity dependence of the isothermal compressibihty, k, defined 
by (H) of the hard-sphere lattice gas on the simple cubic lattice of V 
sites under fully periodic boundary conditions. V = L x L x L; L = 10: 
O, 12: +, 14: □, 16: x, 18: A, 20: 22: o, 24: •. The solid curves 
are obtained by the smoothing procedure of the fourth-order 5-spline. 
Errors are less than the symbol size. 

Figure 2 Activity dependence of the staggered compressibility, defined 
by (p. The meaning of the symbols and the curves is the same as in 
figure 1. Errors are less than the symbol size. 

Figure 3 Activitiy dependence of the fourth-order cumulant, U, defined by 
(1^). The meaning of the symbols and the curves is the same as in 
figure 1. 

Figure 4 Size dependence of the effective critical activities, z^^^{L): +, 
z^.^^{L): o, and ^^^^^(L): x. Errors are less than the symbol size for 
z^.^^{L) and z^^^^^{L). The horizontal line denotes z = 1.0588. 
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Figure 1 
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